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Abstract 

Within the numerically exact solution to the Dicke model proposed previously, we study the 
quantum criticality in terms of the ground-state (GS) energy, fidelity, and the order parameter. 
The finite size scaling analysis for the average fidelity susceptibility (FS) and second derivative of 
GS energy are performed. The correlation length exponent is obtained to be v = 2/3, which is the 
same as that in Lipkin-Meshkov-Glick model obtained previously, suggesting the same universality. 
It is observed that average FS and second derivative of GS energy show similar critical behavior, 
demonstrating the intrinsic relation in the Dicke model. The scaling behavior for the order param- 
eter and the singular part of the GS energy at the critical point are also analyzed and the obtained 
exponents are consistent with the previous scaling hypothesis in 1/N expansion scheme. 

PACS numbers: 42.50.Nn, 64.70.Tg, 03.65.Ud 
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I. INTRODUCTION 



The Dicke model [lj describes the interaction of N two- level atoms (qubits) with a single 
bosonic mode, which is a fundamental model in quantum optics. It exhibits a "superrandi- 
ant" quantum phase transition (QPT) 2J in the thermodynamic limit. In recent years, the 
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Dicke model has attracted considerable attentions 

hand, the quantum entanglement JjJ and Berry phase [l3| have been used to characterize 
the QPTs. On the other hand, the Dicke model is closely related to many recent interesting 
fields in quantum optics and condensed matter physics, such as the superradiant behavior 



of optical cavities 



by an ensemble of quantum dots [ijj] and Bose-Einstein condensates [l5j, and coupled arrays 



The scaling exponents obtained at the critical points are of significance to distinguish 
the universality class of the QPTs among various models. By a modified Holstein-Primakoff 
approach, Vidal and Dusuel has predicted theoretically the nontrivial scaling exponent for 
several quantities in the Dicke model|8l. To our knowledge, the finite-size studies in the Dicke 
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in small size system N < 35, the adiabatic approximation [6J. Recently, by using extended 
bosonic coherent states, a numerically exact technique to solve the Dicke model for large 
system size was developed by the present authors [111]. Therefore, numerical calculations of 
scaling exponents for some key quantities based on convincing treatment for large system 
size are clearly desirable for confirmation. 

Recently, a concept in quantum information theory, i.e. the fidelity has been extensively 
used to identify the QPTs in various man y-b ody systems from the perspective of the ground- 



state (GS) wave functions]! 




271 ] (For more details, please 



refer to a review article [281] ). In a mathematical sense, the fidelity is the overlap between 
two ground states where the transition parameters deviate slightly. However, the fidelity 
depends on a arbitrary small amount of the transition parameters, which in turn yields an 



artificial factor. Zanardi et al 



19[ ] introduced the Riemannian metric tensor and You et al [3] 



proposed the fidelity susceptibility (FS) to avoid this problem independently. The leading 
term of the fidelity are focused to account for the singularity of QPTs in both methods. It 



is thus implied that the fidelity have no singular behavior in t 



transition, which can distinguish different transition types 21 



re Kosterliz-Thouless phase 
22, be|. In addition, the 
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intrinsic relation between the GS fidelity and the derivative of GS energy has been studied 



and it is observed that they play a equivalent role in identifying the QPTs 



221 ] . Since no a 



priori knowledge of the order parameter is needed, it may be a great advantage to study the 
scaling behavior of FS to characterize the universality in quantum critical phenomena [2 ll]. 

In this paper, we extended our previous numerically exact technique to calculate the GS 
energy, FS, and the order parameter in the Dicke model for finite N. The scaling behavior of 
FS, the second derivative of GS energy, and the order parameter of the QPT will be analyzed. 
We will also study the singular part of the GS energy, which is a further correction to the 
regular part sj]. 



II. MODEL HAMILTON! AN 

Without the rotating-wave approximation, the Hamiltonian of N identical atoms inter- 
acting with a single bosonic mode reads 

2A 

H = uja + a + u J z H =(a + + a)J x , (1) 

VN 



where a + and a are the bosonic annihilation and creation operators, ujq and u are the 
transition frequency of the qubit and the frequency of the single bosonic mode, A is the 
coupling constant. J x and J z are the usual angular momentum. There is a conserved parity 
operator II = e l7T ( Jz + N / 2 + a+a ) ) which commutes with the Hamiltonian (1). 



In our previous numerically exact approach the wave function can be expressed in 
terms of the basis {\(p n )b<S) \h n )} where \j,n) is the Dicke state with j = N/2 and \(p n ) b 
is the bosonic state. The latter is given by 

N tr , 

>* = E c ^^rM + + 9n) k e-^ a+ -^ 2 \0) a , (2) 
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where g n = 2\n/u^/N, N tr is the truncated bosonic number in the Fock space of the new 
operator A = a + g n , the coefficient c n ^ can be determined through the exact Lanczos 



diagonalization 291 ] . 



The driving Hamiltonian in the Dicke model is 



Hi = -^={a + + a)J x , (3) 
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and the transition parameter is A. According to the definition 20j, the FS is given by 



where \l/ n is the n-th eigen-states of the Hamiltonian (1). In terms of Eq. (2), we have 

2 TV 

(*„|ffl|¥o> = -^^(y-m)^^.[V7TT^ +1 

m,i,j 

After the FS is well defined above, the fidelity is readily obtained through an artificially 
introduced small amount of the transition parameter 8X 



F(X,5X) = y/1 -5\ 2 S F {\) (5) 



Although the FS might not be a intensive quantity in some cases[2l|], for convenience, we 
will discuss the average FS: xf = Sp(X)/N throughout the paper. It is interesting to note 
that the FS can be obtained through the wave function without the knowledge of the order 
parameter. 

III. FINITE-SIZE SCALING ANALYSIS 

We have studied the leading finite-size corrections to the GS energy in the Dicke 



model llj]. For convenience, we have introduced one dimensionless parameter D = ujq/uj. 
The QPT occurs at the critical point A c = 0.5 in the thermodynamic limit (jj. The scaling 
exponent of the leading finite-size corrections is obtained to be —1.0 ± 0.02 for three typ- 
ical different values D = 0.1, 1, 10, consistent with that by a modified Holstein-Primakoff 



approach [8]. In this paper, we first study the scaling behavior of the singular part of the GS 



energy, which can be obtained by subtracting the leading contribution 8j cO + cl/N, with 
Co = —ujq/2 and c\ = |[— u — ujq + (to 2 + ^o) 1 ^ 2 ] fr° m the total GS energy. We exhibit the 
curve of eo — (co + %) as a function of N at the critical point for different typical values of D 
in Fig. 1. One can obviously observe that the scaling exponent estimated from the slope in 
all cases are —1.325 ±0.005, which is very close to —4/3, confirming the previous theoretical 
prediction in a modified Holstein-Primakoff approach 
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FIG. 1: Scaling of the singular part of the GS energy eo — (co + %>) as a function of N at the critical 
point for D = 0.1, 1, and 5. The inset shows the corresponding slope versus 1/N. 

Next, we illustrate the scaling behavior of the average FS. The finite-size scaling ansatz 
for the average FS to analyze the QPT take the form 27] 



^ = f[N"(\ - A max )] (6) 

Xf 

where x™ ax * s the value of average FS at the maximum point A max , / is the scaling function 
and v is the correlation length critical exponent. This function should be universal for large 
N in the second-order QPTs, which is independent of the order parameter. As shown in 
Fig. 2, with v = 2/3, an excellent collapse in the critical regime is achieved according to 
Eq.(6) in the curve for different large size for three typical values of D. Beyond the critical 
regime, the collapse becomes poor. As N increases, the critical regime become wider. It is 
demonstrated that v is a universal constant and does not depended on the parameter D. It 
is also implied from the collapse that the correlation length behaves like £ oc |A — A c |~ 2 ^ 3 . 
It is very interesting to note that the value of v is the same as that in the Lipkin-Meshkov- 



Glick (LMG) model obtained analytically 30j and numerically [27J], suggesting the same 
universality in the Dicke and LMG model. 

Fig. 3 shows the scaled average FS at the maximum point as a function of N for different 
values of D in log-log scale. A power law behavior x™ ax oc exists in the large N. The 
finite-size exponents extracted from all curves tend to a converging value fx = 0.33 ± 0.02. 
This exponent is also independent of the value of D, and then is a universal constant. 
It is also very close to that in the LMG model, providing another evidence of the same 
universality class of these two models. 

To exhibit the overall properties of the average FS in the whole coupling regime, we 




FIG. 2: Finite-size scaling of the average FS according to Eq. (6) at the critical point for (a) 
D = 0.1, (b) D = 1, and (c) D = 5. 



O 3 



(0=1 








■ - — S--- — 








o D=0.1 ' 






□ D=1.0 






D=5.0 



10 



iog 2 (N) 



13 



FIG. 3: Scaling of the maximum of the average FS as a function of N at the critical point for 
D = 0.1, 1, and 5. The inset shows the corresponding slope versus 1/N. 



calculate the average FS as a function of A for different size. From Fig. 4, we can see 
that S p = Nxf is a intensive quantity when A > A c , while \f is a intensive quantity when 
A < A c . It is pointed out in Ref. 



21 



27| that the intensive average FS in the thermodynamic 
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FIG. 4: The average FS as a function of A at D = 1 for different system sizes. The inset shows the 
collapse of 1 — Xf/x™* as a function of iV 2 / 3 (A — A max ). 

limit scales generally like 

in the vicinity of the critical point. Within the similar analysis, we can get the exponent a = 
1/2 when A > A c and a = 2 when A < A c through the relation a = fi/u, which is readily 
derived from the above scaling ansatz. Also it is interesting to note that the average FS 
becomes divergent as the system size increases, demonstrating a Landau-type transitions in 
the Dicke model. 

To illustrate the intrinsic relation between the ground-state fidelity and the ground state 
energy e$ in the Dicke model, we calculate the second derivative of the GS energy for different 
lattice size at D = 1, which are shown in Fig. 5. Interestingly, similar to the average FS, 
the second derivative of the GS energy also shows divergent behavior around the critical 
point with increasing system size. This divergent behavior is also similar to that observed 



in ID transverse field Ising model 22|. The only difference is that the maximum point X ma x 
approaches the critical point with the increasing system size in the Dicke model more fast 
than in the ID transverse field Ising model, due to the smaller value of the correlation length 
exponent v in the former model. 

Because the scaling A — A max a jV~ 2//3 holds, we try to plot the dimensionless quantities 
1 — i9/'# max , where is either the average FS or the second derivative of the GS energy for 
different size as a function of iV 2//3 (A — A max ) in the insets of Figs. 4 and 5. It is interesting 
to observe that the second derivative of the GS energy exhibits the same scaling behavior 
around the critical point as the average FS. It follows that these two quantities play the 
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FIG. 5: The second derivative of the GS energy as a function of A at D = 1 for different system 
sizes. The inset shows the collapse of 1 — f^/§3pU m as a function of iV 2 / 3 (A — X max ). 

same roles in the characterizing the QPTs in the Dicke model. In the ID transverse field 
Ising model, the value of the critical exponent v = 1, which is equal to that in the classical 
two-dimensional Ising model, the scaling A — A max oc N~ x should be satisfied, which was 
just used in the Figs. 1 and 2 in Ref. |22j ]. 

We then perform the finite size scaling analysis on the order parameter of the QPT, i.e. 
the expectation value of the photon number per atom in the ground-state {afc^/N. In the 
thermodynamic limit, this quantity changes from zero to finite value smoothly when crossing 
the critical point. In Fig. 6, we present this quantity as a function of N for different values 
of D in log- log scale. Derivative of these curves are plotted in the inset. The exponent of the 
order parameter is estimated to be —0.66 ± 0.02, which is consistent well with the predicted 
value 2/3 by using the diagonalizing a expanded Hamiltonian at order 1/N based on 
Holstein-Primakoff representation. 

Finally, to describe the new numerically exact approach with extended coherent states 
in more intuitive detail and see how it connects to results in the thermodynamic limit, 
we calculate the finite but large N wave functions, as done in Ref. jsj]. We calculate the 
ground-state wave function ty(x,y) in the x-y representation for N = 1000 and D — 1 at 
various A which cover A c (iV = 1000). A c (iV = 1000) falls in the range of [0.5070,0.5071], 
according to the position of the peak for FS. After tedious calculations, the numerical results 
for the wave function are shown in Fig. 7, where the displacement is not removed. One can 
observe that the wave packet becomes stretched in a direction with an angle around 7r/4 as 
A increases. This stretching increases up to X C (N), where the wave function is not divergent 
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FIG. 6: Scaling of the order parameter (a^a)/N as a function of iV at the critical point for D = 
0.1, 1, and 5. The inset shows the slope versus 1/N. 

due to the finite-size effect. Exceeding A c , the wave function splits into two peaked one. It 
is also a finite-size effect and not shown in the thermodynamic limit [cf. Fig. 6 in Ref. j^]]. 
With the further increase of A, the two lobes start mixing, a single-peaked wave function is 
then formed like in the thermodynamic limit Q]. It follows that the essential feature in the 
thermodynamic limit already appears for large N, say iV = 1000. Note that in Fig. 12 of 
Ref. [sj], two lobes still separate even at A = 0.7 for j — 5 (i.e. N = 10), different from the 
present observation for iV = 1000. 



IV. CONCLUSION 



In summary, based on our previously proposed numerically exact technique in the finite- 
size the Dicke model, we study the quantum criticality in terms of the GS energy, fidelity, 
and the order parameter. By subtracting the regular part of the GS energy, we perform the 
analysis for the scaling behavior of for the singular part of the GS energy at the critical point. 
The obtained exponent is very close to 4/3, which agree well with the previous theoretical 
prediction in 1 /N expansion scheme based on the Holstein-Primakoff transformation. Then, 
we perform the finite-size scaling analysis for the average FS. Several scaling exponents are 
obtained: v = 2/3, fi = 4/3, and a = 1/2 when A > A c and a = 2 when A < A c . All this 
exponents are the same as those in Lipkin-Meshkov-Glick model obtained both analytically 
and numerically, suggesting the same universality of these two models. We further study the 
scaling behavior of the average FS and the second derivative of GS energy, and observed that 
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FIG. 7: The ground-state wave function ^(x, y) in the x-y representation for N = 1000 and D = 1 
at A = 0, 0.5, 0.507, 0.5071 (from left to right in the upper row), 0.525, 0.55, 0.555, 0.7 (the low row). 
Note that 0.5070 < A c (iV = 1000) < 0.5071. 

these two quantities play the same roles in the QPTs in the Dicke model. Brief comparisons 
with the ID transverse field Ising model are also carried out. We also perform the analysis 
of the finite size scaling effect of the order parameter, and observe that the order parameter 
vanishes as N~ 2 ^ 3 at the critical point, consistent with the previous theoretical prediction. 
Finally, we calculate the finite but large N wave functions. It is observed that the essential 
feature in the thermodynamic limit has already shown up in large system size (N = 1000). 

V. ACKOWNLEDGEMENTS 

The authors acknowledge useful discussion with S.J. Gu and J. Vidal. This work was sup- 
ported by National Natural Science Foundation of China, PCSIRT (Grant No. IRT0754) in 
University in China, National Basic Research Program of China (Grant No. 2009CB929104), 
Zhejiang Provincial Natural Science Foundation under Grant No. Z7080203, and Program 
for Innovative Research Team in Zhejiang Normal University. 



10 



* Corresponding author. Email:qhchen@zju. edu.cn 



[1] R. H. Dicke, Phys. Rev. 93, 99(1954). 

[2] S. Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, England, 
2000). 

[3] C. Emary and T. Brandes, Phys. Rev. E 67, 066203(2003); Phys. Rev. Lett. 90, 044101(2003). 

[4] N. Lambert, C. Emary, and T. Brandes, Phys. Rev. A. 71, 053804(2005). 

[5] N. Lambert, C. Emary, and T. Brandes, Phys. Rev. Lett. 92, 073602(2004). 

[6] G. Liberti, F. Plastina, and F. Piperno, Phys. Rev. A 74, 022324 (2006). 

[7] V. Buzek, M. Orszag, and M. Rosko, Phys. Rev. Lett. 94, 163601(2005). 

[8] J. Vidal and S. Dusuel, Europhys. Lett. 74, 817(2006). 

[9] J. Reslen, L. Quiroga, and N. F. Johnson, Europhys. Lett. 69, 8(2005). 

[10] F. Plastina, G. Liberti, and A. Carollo, Europhys. Lett. 76, 182(2006). 

[11] Q. H. Chen, Y. Y. Zhang, T. Liu, and K. L. Wang, Phys. Rev. A 78, 051801(R) (2008). 

[12] A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature (London) 416, 608(2002); T. J. Osborne 

and M. A. Nielsen, Phys. Rev. A. 66, 032110(2002)). 

[13] A. C. M. Carollo and J. K. Pachos, Phys. Rev. Lett. 95, 157203(2005); S.- L. Zhu, Phys. 

Rev. Lett. 96, 077206(2006). 

[14] M. Scheibner et al., Nature Phys. 3, 106(2007). 

[15] D. Schneble et al., Science 300, 475 (2003). 

[16] M. J. Hartmann et al., Nature Phys. 2, 849(2006); A. D. Greentree et al., ibid. 2, 856(2006). 

[17] H. T. Quan, Z. Song, X. F. Liu, P. Zanardi, and C. P. Sun, Phys. Rev. Lett. 96, 140604 
(2006). 

[18] P. Zanardi, and N. Paunkovic, Phys. Rev. E. 74, 0331123 (2006). 

[19] P. Zanardi, P. Giorda, and M. Cozzini, Phys. Rev. Lett. 99, 100603(2007). 

[20] W. L. You, Y. W. Li, and S. J. Gu, Phys. Rev. E. 76, 022101 (2007). 

[21] S. J. Gu, H. M. Kwok, W. Q. Ning, and H. Q. Lin, Phys. Rev. B. 77, 245109(2008). 

[22] S. Chen, L. Wang, Y. J. Hao, and Y. P. Wang, Phys. Rev. A. 77, 032111 (2008). 

[23] H. Q. Zhou, R. Orus, G. Vidal, Phys. Rev. Lett. 100, 080601 (2008); H. Q. Zhou, J. P. 



Barjaktarevic, |arXiv:cond-mat/0 701608 H. Q. Zhou, J. H. Zhao, B. Li. larXiv:0704~2 940 , 

11 



[24] J. Ma, L. Xu, H. N. Xiong, and X. G. Wang, Phys. Rev. E. 78, 051126 (2008). 
[25] K. W. Sun, Y. Y. Zhang, and Q. H. Chen, Phys. Rev. B. 79, 104429 (2009). 
[26] L. Campos Venuti, M. Cozzini, P. Buonsante, F. Massel, N. Bray-Ali, P. Zanardi, Phys. Rev. 
B 78, 115410 (2008) 

[27] H. M. Kwok, W. Q. Ning, S. J. Gu, and H. Q. Lin, Phys. Rev. E. 78, 032103 (2008). 
[28] S. J. Gu, arxiv:quant-ph/0811.3127 

[29] G. S. Tian, L. H. Tang, and Q. H. Chen, Europhys. Lett. 50, 361(2000); Phys. Rev. B63, 
054511(2001). 

[30] S. Dusuel and J. Vidal Phys. Rev. Lett. 93, 237204 (2004); Phys. Rev. B 71, 224420 (2005). 



12 



